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ABSTRACT 

Prediction of the Sun's magnetic activity is important because of its effect on space environment 
and climate. However, recent efforts to predict the amplitude of the solar cycle have resulted in 
diverging forecasts with no consensus. Yeates et al. (2008) have shown that the dynamical memory 
of the solar dynamo mechanism governs predictability and this memory is different for advection- 
and diffusion-dominated solar convection zones. By utilizing stochastically forced, kinematic dynamo 
simulations, we demonstrate that the inclusion of downward turbulent pumping of magnetic flux 
reduces the memory of both advection- and diffusion-dominated solar dynamos to only one cycle; 
stronger pumping degrades this memory further. Thus, our results reconcile the diverging dynamo- 
model-based forecasts for the amplitude of solar cycle 24. We conclude that reliable predictions for 
the maximum of solar activity can be made only at the preceding minimum-allowing about 5 years 
of advance planning for space weather. For more accurate predictions, sequential data assimilation 
would be necessary in forecasting models to account for the Sun's short memory. 
Subject headings: magnetic fields - Sun: activity - Sun: dynamo - Sun: interior 



1. INTRODUCTION 

Sunspots are strongly magnetized regions on the so- 
lar surface and have been systematically observed over 
the last four centuries. These observations show that 
the magnetic output of the Sun is roughly periodic, with 
the number of sunspots on the solar surface increasing 
and decreasing in a cyclic fashion with a mean period 
of 11 years. Sunspots are the seats of solar magnetic 
storms-solar flares and coronal mass ejections-that eject 
large amounts of magnetized plasma into space generat- 
ing severe space weather. When these storms are Earth 
directed, they pose a serious hazard to satellites, air- 
traffic on polar routes and modern-day technologies such 
as telecommunication, GPS networks and electric power 
grids at high latitudes (US National Academy of Sciences 
2008). The slow, long-term evolution of magnetic activ- 
ity also affects the radiative output of the Sun (Lean 
1997; Krivova et al. 2007)-the primary natural energy 
input into the climate system, and is therefore important 
from the perspective of climate dynamics. This realiza- 
tion has driven many recent efforts to predict the activity 
of the Sun, however, with disappointingly, diverging re- 
sults (Pesnell 2008). 

The magnetic activity of the Sun is generated through 
a magnetohydrodynamic (MHD) dynamo mechanism in- 
volving interactions between plasma flows and magnetic 
fields within the solar convection zone (SCZ). The cur- 
rent understanding (Charbonncau 2010; Nandy 2012) of 
the solar cycle is that the strong toroidal (</>) compo- 
nent of the solar magnetic field is generated in the solar 
interior through stretching of the poloidal (r-9) compo- 
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nent of the magnetic field by the differentially rotating 
solar plasma. These toroidal flux tubes are magneti- 
cally buoyant and can rise from the interior and erupt 
through the solar surface forming bipolar sunspot pairs. 
Traditionally two mechanisms-namely the mean-field a- 
effect driven by helical turbulence, and the Babcock- 
Leighton process (Babcock 1961; Leighton 1969; Munoz- 
Jaramillo et al. 2010) related to the decay and disper- 
sal of tilted bipolar sunspot pairs at near-surface layers- 
are thought to recreate the poloidal component back 
from the toroidal field. Here, we focus on the Babcock- 
Leighton solar dynamo model as the Babcock-Lcighton 
process can be constrained by surface magnetic field ob- 
servations and long-term observations spanning multiple 
solar cycles suggests that the Babcock-Leighton mecha- 
nism predominantly contributes to the solar poloidal field 
(Dasi-Espuig et al. 2010). 

Since the two components of the Sun's magnetic field 
are produced at spatially segregated layers in the SCZ, 
transport of magnetic flux between them is necessary for 
the functioning of the cycle. The finite time required 
for magnetic flux transport introduces a memory into 
the system which allows for dynamo model based solar 
cycle predictions. While magnetic buoyancy plays the 
dominant role in transporting toroidal field to the sur- 
face, there are diverse competing mechanisms for trans- 
porting the poloidal field generated at near-surface layers 
back to the solar interior. In advection-dominatcd mod- 
els meridional circulation is assumed to be the primary 
mechanism, while in diffusion-dominated models turbu- 
lent diffusion is assumed to be the primary mechanism 
for achieving the transport of poloidal field into the gen- 
erating layer for the toroidal component. Amongst the 
many predictions (Pesnell 2008) till date for the ampli- 
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tude of solar cycle 24, the only two physics (dynamo) 
based forecasts disagree strongly. One predicts a very 
strong cycle (Dikpati et al. 2006) while the other a very 
weak cycle (Choudhuri et al. 2007). Yeatcs et al. (2008) 
have shown that the diverging predictions result from 
different assumption by these predictive models regard- 
ing the nature of the SCZ; while Dikpati et al. (2006) 
assume a advection-dominatcd SCZ, Choudhuri et al. 

(2007) assume a diffusion-dominated SCZ. Yeates et al. 

(2008) have demonstrated that the memory of the solar 
dynamo is longer in the advection-dominated case (last- 
ing over several solar cycles) as opposed to the diffusion- 
dominated case (lasting over one cycle) and this they 
argue is the origin of the disagreement in the predictions 
(see also Jiang et al. 2007). Solar observations can- 
not conclusively determine whether the Sun's convection 
zone is diffusion-dominated or advection-dominatcd, fur- 
ther compounding the problem of solar cycle predictions. 

However, another important flux transport mechanism 
in the SCZ, namely the turbulent pumping of magnetic 
flux, has not received adequate attention and has been 
completely ignored in solar cycle prediction models. The- 
oretical analysis and magneto-convection simulations by 
several authors (Petrovay & Szakaly 1993; Brandenburg 
et al. 1996; Tobias et al. 1998, 2001; Dorch & Nordlund 
2001; Ossendrijver et al. 2002; Zicgler & Rudiger 2003; 
Kapyla et al. 2006a; Rogachcvskii et al. 2011; Pipin & 
Seehafer 2009; Sechafer & Pipin 2009) suggest that mag- 
netic fields in the SCZ are pumped preferentially along 
radial direction with maximum pumping speeds of a few 
m/s. This radial pumping is particularly effective on 
the weak fields which are not subject to strong buoy- 
ancy forces, mainly acts vertically downward and pre- 
sumably results from a combination of factors including 
the presence of strong, coherent downward plumes and 
radial gradient in turbulence. An equatorward directed 
latitudinal pumping, mainly concentrated near low lati- 
tudes, also arises when the rotation becomes important 
(Ossendrijver et al. 2002; Kapyla et al. 2006a). Many 
studies indicate that turbulent pumping plays a dynam- 
ically important role in the context of the solar and the 
stellar dynamo mechanism (Brandenburg et al. 1992; 
Kapyla et al. 2006b; Guerrero & de Gouveia Dal Pino 
2008; Kitchatinov & Rudiger 2012; Do Cao & Brun 2012; 
Cameron et al. 2012). Nevertheless, the effect of turbu- 
lent pumping on solar cycle memory and predictability 
has been hitherto unexplored. 



2. MODEL 

Here, within the domain of a stochastically forced 
kinematic solar dynamo model, in addition to meridional 
circulation and turbulent diffusion mediated transport 
of the poloidal field, we explore the effects of vertically 
downward turbulent pumping of poloidal magnetic field 
as a dominant flux transport mechanism. The evolution 
of the poloidal and toroidal components of the solar 
magnetic fields in the framework of the (a-51) kinematic 
solar dynamo model are given by 
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with s = rsm9. Here A is the vector potential for the 
poloidal component of the magnetic field (B p =V x A) 
and B is the toroidal magnetic field. The internal rota- 
tion of the Sun denoted by 57 is determined from an an- 
alytic fit to helioseismic data. Meridional circulation is 
given by v = v r e r + veeg; it is directed poleward near the 
surface layers, equatorward near the base of the SCZ and 
is subject to mass conservation (i.e., satisfies V.pv = 0). 
The poloidal field source-term a-which is concentrated at 
near-surface layers-models the Babcock-Lcighton mech- 
anism of poloidal field creation. Effective magnetic dif- 
fusivities for the poloidal and toroidal component of the 
field are denoted by r\ v and r\ t , respectively. We solve the 
above dynamo equations with appropriate choices of pa- 
rameters and boundary conditions within a r-9 domain 
that includes the SCZ and one solar hemisphere. The 
typical set-up of this model has been well-studied (Nandy 
& Choudhuri 2002; Chatterjee et al. 2004; Yeates et al. 
2008) and for ease of comparison-spccifically in the con- 
text of solar cycle memory, here we use the same parame- 
ters as utilized in the study of Yeates et al. (2008)-where 
all input parameters have been described in details. 

Turbulent pumping appears as an advective term in 
the mean field induction equation, is shown to be more 
effective on the weaker poloidal fields, and previous stud- 
ies indicate that it is directed downwards across all lat- 
itudes (Ossendrijver et al. 2002; Kapyla et al. 2006a). 
The stronger toroidal component of the field which exists 
as localized and magnetically buoyant flux tubes in the 
SCZ is expected to be affected less by the pumping mech- 
anism. The memory of the solar cycle-exploring which is 
the primary aim of this study-is also known to be affected 
primarily by the time delay introduced due to the trans- 
port of poloidal flux downward into the solar interior. 
Motivated by these considerations, and to keep our sim- 
ulations and the interpretation of the results tractable, 
we include only the effects of downward turbulent pump- 
ing on the poloidal component of the magnetic field. To 
achieve this we add an extra term in the advective part 
of Eq. 1, given by v p = 7 r e r , to model radial pump- 
ing. Note that this parameterization in not intended to 
model a global downward movement of plasma, but only 
that of magnetic flux and hence does not need to sat- 
isfy mass conservation. Guided by numerical simulations 
(Ossendrijver et al. 2002; Kapyla et al. 2006a) we use 
the following profile for j r : 
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The value of /707- determines the amplitude of "f r , which 
we vary in our simulations to explore its effects on so- 
lar cycle memory. The radial pumping profile in the 
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SCZ based on Eq. 3 is depicted in Fig. 1. It is nega- 
tive (i.e., vertically downward) throughout the convec- 
tion zone and vanishes below 0.7 R (near the base of the 
SCZ). Its value is maximum near the poles at a radius 
around 0.9R and decreases towards the equator and in 
the solar interior. 
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Fig. 1. — Dotted, dashed, dot-dashed and solid lines show the 
variation of radial pumping j r at co- latitudes 6 = 0°, 45°, 60° and 
90° , respectively, with fractional solar radius from the solar surface 
to the solar interior. 



3. SIMULATIONS 

Fluctuations in the solar cycle in the context of the 
Babcock-Leighton dynamo arise from the random buf- 
feting of rising magnetic flux tubes, resulting in a scat- 
ter around the mean tilt angle of bipolar sunspot pairs 
given by Joy's law-which is in fact observed (Dasi-Espuig 
et al. 2010). Thus, to imbibe the stochastic nature of 
the Babcock-Leighton poloidal field generation process 
in our model, we introduce random fluctuations in the en- 
effect by setting a = a hllsc + afi uc a(t,T COI ), with 

Cubase — 

30 m s . For our simulations we take afluc = Cubase (i-e., 
with fluctuations up to 100%). Here, a is an uniform ran- 
dom deviate between [—1 and +1]. The coherence time 
r cor is chosen such that there are around ten fluctuations 
in a (sunspot) cycle. Treating stochastic fluctuations in 
this manner is standard practice, even in dynamo mod- 
els based on the mean-field a-effect models (Pipin et al. 
2012) . Although a precise observational constraint on the 
amplitude of the Babcock-Leighton a-effect docs not ex- 
ist, an order of magnitude estimate is possible based on 
the efficacy of the decay and dispersal of active region 
flux governed by diffusion and meridional circulation. 
For the diffusive dispersal velocity (v v = r]/L), consider- 
ing an active region scale L = 10 7 m and rj = 1 x 10 8 
m 2 s _1 , yields v ri = 10 m s . In conjunction, a typical 
near-surface meridional circulation speed on the order of 
20 m s^ 1 motivates our choice of the Babcock-Leighton 
Qf-cffcct amplitude. 

With stochastic forcing, wc explore the effect of turbu- 
lent pumping on solar cycle memory through simulations 
in two different regimes, diffusion-dominated (vo = 15 
m s" 1 , t/o = 1 x 10 12 cm 2 s" 1 ) and advection-dominated 
regime (vq — 26 m s" 1 , 170 = 1 x 10 12 cm 2 s _1 ). With 
these parameter choices, in the former case turbulent dif- 
fusivity plays a dominant role, while in the latter case ad- 
vection by meridional circulation plays a dominant role 
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Fig. 2. — Representative butterfly diagram of the solar toroidal 
field at the base of the convection zone (solid and dashed contours 
depict positive and negative fields, respectively) is overlaid on the 
radial field at the solar surface (in grayscale, with white denoting 
positive and black denoting negative values). This diagram depicts 
the spatio-temporal field evolution from a stochastically forced dy- 
namo simulation in the advection-dominated regime with a pump- 
ing speed of 2 m s —1 . These simulations exhibit cyclic reversals as 
well as fluctuations in the magnetic field amplitude. 

in the transport of poloidal flux from the upper to the 
lower layers of the SCZ (Yeates et al. 2008). We vary 
the maximum pumping speed progressively from to 4 
m s _1 in both these regimes, allowing transport by tur- 
bulent pumping to become increasingly important in the 
dynamics of magnetic flux transport. 

4. RESULTS 

Figure [2] shows a butterfly diagram of the toroidal field 
at the base of the convection zone overlaid on the surface 
radial field from a simulation in the advection-dominated 
regime, with a pumping speed of 2 m s _1 . The ampli- 
tude of the generated toroidal and poloidal fields vary 
from one cycle to another in these stochastically forced 
runs. Over continuous simulations spanning multiple so- 
lar cycles we calculate the polar radial flux <\> r located 
within the latitude 70° to 89° at the solar surface and the 
toroidal flux <p tor located within a layer r = O.677i?0- 
O.726i? and latitudinal extent 10°-45° - i.e., in the 
low latitude tachocline. From this we calculate the peak 
value of 4>tor and 4> r of each cycle to compute the cor- 
relation between the peak surface radial flux for cycle n 
with the toroidal flux of subsequent cycles. The extent 
of this correlation, as argued by Yeates et al. (2008), es- 
tablishes the dynamical memory of the solar cycle. Note 
that this approach does not have to rely on identifying 
the precise epoch of solar maximum or minimum based 
on any qualitative appraisal of the butterfly diagram and 
is therefore both quantitative and objective. The simu- 
lated toroidal flux for a cycle is on the order of 10 25 
Mx which is close to what is observed but the simulated 
peak polar flux of about 10 24 Mx is higher than the ob- 
served polar flux at solar minimum of about 10 22 Mx. 
However, the latter observations are based on relatively 
low-resolution magnetograms and the recent high resolu- 
tion observations of unipolar kilogauss flux tubes in the 
polar region (Tsuneta et al. 2008) can potentially rec- 
oncile this (for detailed arguments see Choudhuri 2003 
and the online supplementary information in Nandy et 
al. 2011). 

The influence of increasingly dominant turbulent 
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TABLE 1 

Correlation coefficients (r 3 ) and percentage significance 

LEVELS (p) FOR PEAK SURFACE RADIAL FLUX <3> r OF CYCLE n 
VERSUS PEAK TOROIDAL FLUX <3>tor OF DIFFERENT CYCLES FOR 275 
SOLAR CYCLES DATA. THE FIRST COLUMN DENOTES THE 
AMPLITUDE OF THE TURBULENT PUMPING SPEED IN VARIOUS 
SIMULATION STUDIES. THE TOP ROW CORRESPONDS TO THE CASE 
WITHOUT TURBULENT PUMPING AND SUBSEQUENT ROWS 
CORRESPONDS TO SIMULATIONS WITH INCREASING PUMPING SPEEDS. 















Uii. Dom. 


Adv. 


Dom. 


Pumping 




Parameters 




r s (p) 




r s ip) 




*r 


("■) 


& *tor 


(n) 




0.19 (99.9) 


0.57 


(99.9) 




*r 


J 1 ) 


& *tor 


>4 


1) 


0.64 (99.9) 


0.77 


(99.9) 


m s — 


*i 


n) 


& *tor 


[n4 


2) 


0.04 (55.9) 


0.46 


(99.9) 




ih 
?r 


n) 


& *tor 


[n4 


" 3) 


n oo /no n\ 
0.Z2 (99. 9) 


0.27 


(99.9) 




*r 


( rt ) 


& *tor 


(n) 




-0.06 (67.0) 


0.41 


(99.9) 




*, 


ji) 


& $tor 


[n4 




0.67 (99.9) 


0.72 


(99.9) 


1 m s -1 


<£>■■ 


n) 


& *tor 


[n4 


2) 


0.09 (83.9) 


0.29 


(99.9) 




*r 


a ) 


& *tor 


[n4 


-3) 


-0.02 (26.5) 


-0.01 


(18.9) 




*r 


(") 


& <J>to r 


(n) 




0.12 (94.9) 


0.19 


(99.8) 




*i 


") 


& *tor 


>4 


1) 


0.43 (99.9) 


0.75 


(99.9) 


2 m s" 1 


4> r 


") 


& *tor 


;n4 


2) 


-0.16 (99.9) 


0.07 


(73.8) 




<3>r 


» ) 


& <I>tor 


[n4 


-3) 


-0.02 (20.8) 


-0.10 


(89.8) 




*r 


(") 


& *tor 


(n) 




0.11 (49.2) 


0.29 


(92.0) 




*r 


ji) 


& $tor 


[n4 


1) 


0.32 (99.9) 


0.62 


(99.9) 


3 m s -1 


*, 


ji) 


& $tor 


[n4 


2) 


-0.18 (99.6) 


0.07 


(78.0) 




*r 


Ji) 


& *tor 


[n4 


'3) 


0.03 (36.6) 


-0.10 


(91.6) 




#r 


[n] 


& $tor 


(n) 




0.19 (99.8) 


0.30 


(99.9) 




*i 


Jl) 


& $tor 


[n4 


1) 


0.26 (99.9) 


0.46 


(99.9) 


4 m s _1 


*i 


Ji) 


& *tor 


[n4 


2) 


-0.16 (99.3) 


0.07 


(72.8) 




*r 


n) 


& *tor 


[n4 


-3) 


-0.10 (91.9) 


-0.22 


(99.9) 



pumping on the memory of the solar cycle is summarized 
in Table 1 for both advection- and diffusion-dominated 
convection zones. Spearman's rank correlation coeffi- 
cients and significance levels have been computed be- 
tween the peak polar and toroidal fluxes of various cycles. 
As a consistency check, for the case without any turbu- 
lent pumping, we first recover the result of Yeates et al. 
(2008) that the memory of the solar dynamo is short 
in the diffusion-dominated regime-with the correlation 
lasting only between the polar flux of cycle (n) and the 
toroidal flux of cycle (n4-l), but long in the advection- 
dominatcd regime-where the poloidal flux of cycle (n) 
correlates with the toroidal flux of cycles (n), (n4T), 
(n4-2) and (n4-3) with the memory degrading slowly. 

However, it is observed that when turbulent pump- 
ing becomes increasingly important, the advcction- 
dominatcd dynamo system starts losing its memory and 
at a modest pumping speed of 2 m s , the signifi- 
cant correlation is only between the polar flux of cy- 
cle (n) and the toroidal flux of the next cycle (n+1). 
Unexpectedly therefore, the diffusion-dominated and 
advection-dominated dynamo regimes become similar- 
at least in terms of memory-when flux transport by 
turbulent pumping becomes significant. Figure [3] (for 
the advection-dominated regime) and Fig. @] (for the 
diffusion-dominated regime) , where we plot the relation- 
ship between poloidal and toroidal flux of various cycles 
with a pumping speed of 2 m s _1 , reinforce this con- 
clusion. Moreover, as summarized in Table 1, we find 
that when the amplitude of turbulent flux pumping is 
made stronger the memory of the solar cycle reduces for 
all cases, and for a pumping speed of 4 m s" 1 even the 
one cycle memory is severely degraded. To check the 



parameter dependence of the results we have performed 
some simulations with changes in meridional circulation 
in conjunction with fluctuations in the poloidal field gen- 
eration process. Preliminary findings indicate that the 
above results remain qualitatively unchanged. 

5. IMPLICATIONS FOR SOLAR CYCLE PREDICTIONS 

The overreaching conclusion from this study is that the 
memory of the solar cycle is short, extending over a time 
period barely spanning one solar cycle (between the mini- 
mum in activity of one cycle and the maximum in activity 
of the next). Our analysis, which leads to that conclusion 
is remarkable in that it resolves the ambiguity that ex- 
isted between solar cycle predictions based on the differ- 
ing memory of advection- and diffusion-dominated solar 
dynamos. Expressed succinctly, the relative efficiency 
of advective or diffusive flux transport, by meridional 
circulation or turbulent diffusion, respectively, becomes 
less significant in the presence of turbulent pumping of 
magnetic flux and the latter now governs the dynamical 
memory of the system. This is because turbulent pump- 
ing dominates the vertical flux transport by coupling the 
two spatially segregated source layers of the solar dy- 
namo. 

A discussion on the timescales corresponding to com- 
peting flux transport processes in the SCZ illuminates 
the physics underlying our results. A meridional flow 
speed of vq = 26 m s results in a advective flux trans- 
port timescale (L/vq, with L estimated along a typical 
stream-line looping through the pole) ranging from fO 
to 20 yr to transport the surface poloidal field to the 
mid-to-low latitude tachocline-where the toroidal field 
spawns sunspots of the next cycle. Turbulent diffusion 
with r)o = Ix 1 12 cm 2 s _1 results in a diffusion timescale 
(L 2 /?7o, where L signifies the depth of the SCZ) of 14 yr. 
On the other hand, even a modest turbulent pumping 
velocity of j r — 2 m s _1 results in a comparably much 
shorter, pumping timescale (L/j r , where L signifies the 
depth of the SCZ) of about 3.4 yr. Thus, the dynamical 
memory of the solar cycle becomes very short on inclu- 
sion of turbulent pumping. Faster pumping results in 
an even shorter memory, which becomes comparable to 
the coherence time of stochastic fluctuations in the near- 
surface poloidal field generation process. This results in 
fluctuations of the poloidal field, say between the mini- 
mum of cycle (n) and the maximum of cycle (n4-l), be- 
ing reflected in the tachocline where the toroidal field of 
cycle (n4T) is stored and amplified, thereby, degrading 
even the one cycle memory of the system. 

In summary, we find that the solar cycle memory is 
short, extending only over one cycle; the poloidal field 
at the minimum of a given cycle may be used to predict 
the amplitude of the next sunspot cycle. This allows for 
space weather preparedness and planning about 5 years 
in advance, but presumably no longer-at least based on 
this current understanding. Nonetheless, even one cy- 
cle forecasts may be fraught with difficulties as efficient 
turbulent flux pumping from the surface to the solar in- 
terior can degrade this one cycle memory. Perhaps, in 
retrospect, this explains why multiple attempts-even the 
precursor-based ones-to forecast the amplitude of solar 
cycle 24 have diverged. In the light of these findings, 
we believe that the best strategy for accurate solar cy- 
cle forecasts lies in data assimilation techniques that can 
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Fig. 3. — Relationship between the peak (polar) radial flux <j> T {n) and the peak toroidal flux iptor of cycle (a) n (b) n + 1, (c) n + 2, and 
(d) ri + 3 in the advection-dominatcd regime with a pumping speed amplitude of 2 m s -1 . The flux values are in units of 10 25 Mx. The 
Spearman's rank correlation coefficients (r s ) along with significance levels are inscribed. The dynamical memory persists only between the 
peak polar flux of cycle (n) and the peak toroidal flux of cycle (n+1), i.e., spanning a time which is less than one solar cycle. 
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compensate for the loss of memory through data assim- 
ilation at frequent intervals in forecasting models, with 
the requirement that those intervals be relatively shorter 
than the memory of the solar cycle. 
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